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INPUT SYSTEM BASED ON A THREE-DIMENSIONAL INERTIAL NAVIGATION 
SYSTEM AND TRAJECTORY ESTIMATION METHOD THEREOF 

CROSS-REFERENCE TO RELATED APPLICATIONS 
This application claims the priority of Korean Patent Application No. 2003-16779, 
filed on March 13, 2003, in the Korean Intellectual Property Office, the disclosure of which is 
incorporated herein in its entirety by reference. 

BACKGROUND OF THE INVENTION 

1. Field of the Invention 

An apparatus consistent with the present invention relates to an input system based on 
a three-dimensional inertial navigation system and a trajectory estimation method thereof. 

2. Description of the Related Art 

With progress in the information age, it is more frequently demanded to access and 
produce information anytime and anywhere, rather than a specific space in computers- 
equipped environments. Such demands are satisfied with personal portable devices that have 
been developed and are widely spread, but further research is still required on input systems 
suitable for such personal portable devices. The specifications required for the input system 
for personal portable devices include an easy-to-carry and easy-to-input device, thus, there is 
a demand to develop an input system that is smaller than personal portable devices as well as 
being easy-to-carry. Further, the nature of the personal portable devices is to facilitate the 
inputs of new data anytime and anywhere, so that an input system is required which enables 
users to input such data as they naturally write. In particular, if such an input system can 



recover writing motions as naturally done on a general plane, a free space, or a sheet of paper 
so as to enable users to input characters, graphics, gestures, and so on, the input system can 
be utilized as an excellent input system since it is excellent in the aspect of "universal 
purpose" and also does not require a particular learning course for use. 

[04] Accordingly, in order to meet the demands on the input system as described above, an 

input system has been proposed that can input data in space using the concept of a three- 
dimensional inertial navigation system. 

[05] The three-dimensional inertial navigation system refers to a system that detects three- 

axis acceleration information and three-axis angular velocity information and applies the 
detected information to a moving object in order to estimate positions and attitudes of the 
moving object in the three-dimensional space. The three-dimensional inertial navigation 
system calculates the attitudes of the object and corrects the acceleration information using an 
integral value of angular velocity information estimation values, accordingly, calculates 
velocity information by integrating the corrected acceleration information once, and 
calculates position information by integrating the corrected acceleration information twice. 

[06] Fig. 1 is a block diagram illustrating an input system using a conventional inertial 

navigation system. 

[07] The input system using the conventional inertial navigation system includes an input 

part 10 and a host device 20. 
[08] The host device 20 displays an image corresponding to motions based on the motions 

from the input part 10. 

[09] The input part 10 includes acceleration sensors 11, angular velocity sensors 12, a 

rotation angle information-computing unit 13, conversion-computing unit 14, and a 
transmission unit 15. 
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[10] The acceleration sensors 1 1 generate acceleration information A b [A bx , A by , A bz ) 

corresponding to motions based on the motions of the input part 10, and output the same to 
the conversion-computing unit 14. The subscript b in the acceleration information denotes 
acceleration information for a body frame of the input part and the subscripts x, y 9 and z 
denote information on x-axis, y-axis, and z-axis. 

[11] The angular velocity sensors 12 generate angular velocity information 

W btybx> W by> W bz ) corresponding to motions based on the motions of the input part 10, and 
output the same to the rotation angle information-computing unit 13. The subscript b of the 
angular velocity information denotes angular velocity information for the body frame of the 
input part, and the subscripts Jt, y, and z denote information on x-axis, y-axis, and z-axis. 

[12] The rotation angle information-computing unit 13 inputs angular velocity information 

from the angular velocity sensors 12. The inputted angular velocity information is converted 
into rotation angle information 0, \|/) through a predetermined computing process. The 
predetermined computing process for the conversion into the rotation angle information is 
well known to those skilled in the art, so a detailed description on the process will be omitted 
in the disclosure. 

[13] The con version -computing unit 14 inputs acceleration information and rotation angle 

information. The conversion-computing unit 14 calculates the attitudes of the input part 10 
based on the rotation angle information, corrects the acceleration information based on the 
rotation angle information, calculates velocity information by integrating the corrected 
acceleration information once, and calculates position information by integrating the 
corrected acceleration information twice. Such calculations can be performed based on 
Equations 1, 2, and 3 as follows. 
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[14] 



Prior to explanations of the following Equations, definitions are made for terms used 



in the Equations as follows: 



[15] 



• A body frame is a frame from which acceleration information and angular velocity 



information can be detected in correspondence to actual motions of the input part 10. 



[16] 



• A navigation frame is a reference frame for calculating information applicable for a 



host device 20 by applying predetermined computing matrices based on a result detected 
from the body frame. Herein, the coordinate axes for the navigation frame are referred to as 
X 0 , Y 0 , and Zo. 



frame and the body frame coincide with each other when there is no rotation at an initial time. 
When the navigation frame and the body frame coincide with each other at the initial time 
and the coordinate axes for the body frame are referred to as X 0 , Y 0 , and Zo, VP denotes the 
rotation angle information for the Zo axis. Further, 9 denotes rotation angle information for a 
Yi axis indicating an axis after the Y 0 axis is rotated as much as *F. <D denotes rotation angle 
information for a X 2 axis indicating an axis after the X 0 axis is rotated as much as *F and G, 
respectively. 

[18] • A 3 x 3 matrix for converting a vector in the body frame into a vector in the 

navigation frame is as follows: 



[17] 



• Rotation angle information is as follows: The coordinate axes for the navigation 



[19] 




respectively, and the % denotes rotation angle information in a 3 x 1 matrix. 



[20] 



• The following is a matrix for calculating the rotation angle information based on the 



angular velocity information in the body frame: 
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R(Z)= 0 </> c -<t> s 

_o 4>je e <p c ie c \ 

Here, the subscripts s, c, and t stand for sine, cosine, and tangent functions, 
respectively, and the^ denotes the rotation angle information in a 3 x 1 matrix. 

If the A b denotes the acceleration information in the body frame, the W b denotes the 
angular velocity information in the body frame, R denotes a matrix for calculating the 
rotation angle information based on the angular velocity information in the body frame, G is 

the gravitational acceleration, and the Q n b denotes a matrix for converting a vector in the 
body frame into a vector in the navigation frame, position information in the navigation 
frame is P n , velocity information in the navigation frame is V n , and the rotation angle 

information is % . These variables can be expressed in Equations 1, 2, and 3 as follows: 

[Equation 1] 

dP =V 

n n 

[Equation 2] 

dv n =clhc)\-G 

[Equation 3] 

dz = R(Z)W b 

Here, P n , V n , A b , G, and W b are a 3 x 1 matrix, respectively, G is a matrix of 
[0,0, > and dP n , dV n , and &x denotes differentiations of P n , V n , and x with respect to 
time, respectively. 

Accordingly, the acceleration information and the angular velocity information for the 
body frame outputted from the acceleration sensors and the angular velocity sensors 



respectively are computed and converted into velocity information and position information 
for the navigation frame, and a motion image corresponding to motions in the three- 
dimensional space based on the conversion-computed position information is projected to an 
optimal projection plane, and then outputted to the transmission unit 15. 



[25] The transmission unit 15 outputs the projected position information to the host device 

20. 

[26] By the above fashion, the motions of the input part 10 are displayed on the host 

device 20. 

[27] However, with both the acceleration sensors and the angular velocity sensors used as 



above, the input part becomes heavy, which brings out a drawback that it is not suitable for 
portability. Also, there is a drawback in that the price of the angular velocity sensors is so 
high that the input parts using them become expensive. Further, there exists a drawback in 
that an amount of power consumption becomes large since the angular velocity sensors and 
the acceleration sensors have to be all driven. Furthermore, there exists troublesomeness in 
that a correction job has to be performed in case that the input part is equipped with the 
angular velocity sensors since the use of the angular velocity sensors inevitably requires an 
initial correction job. 

SUMMARY OF THE INVENTION 

[28] Illustrative, non-limiting embodiments of the present invention overcome the above 

disadvantages and other disadvantages not described above. Also, the present invention is 
not required to overcome the disadvantages described above, and an illustrative, non-limiting 
embodiment of the present invention may not overcome any of the problems discussed above. 

[29] An apparatus consistent with the present invention provides an input system based on 

a three-dimensional inertial navigation system and a trajectory estimation method thereof, 

6 



capable of recovering motions in the three-dimensional space by using only three acceleration 
sensors without any angular velocity sensor, that is, without using both three acceleration 
sensors and three angular velocity sensors as in the three-dimensional inertial navigation 
system. 

The above aspect is achieved by providing an input system based on a three- 
dimensional inertial navigation system and having an input part and a host device, which 
detects motion position information corresponding to three-dimensional motions of the input 
part and outputs the detected motion position information to the host device. The input 
system comprises acceleration sensors for outputting acceleration information just before the 
motions (hereinafter, referred to as pre-motion acceleration information where applicable), 
acceleration information on the motions (hereinafter, referred to as motion acceleration 
information where applicable), and acceleration information just after the motions 
(hereinafter, referred to as post-motion acceleration information where applicable); a rotation 
angle information estimation-computing portion for estimating rotation angle information d>, 
G, and *F on the motions (hereinafter, referred to as motion rotation angle information where 
applicable) through a predetermined computing process based on the outputted pre-motion 
acceleration information and post-motion acceleration information; a conversion-computing 
unit for calculating the motion position information based on the estimated motion rotation 
angle information and the outputted motion acceleration information; and an optimal plane- 
computing unit for projecting the motion position information onto an optimal plane. 

The rotation angle information estimation-computing portion includes a first 
computing unit for calculating rotation angle information Ol, 01, and ¥1 just before the 
motions (hereinafter, referred to as pre-motion rotation angle information where applicable) 
and rotation angle information <D2, 92, and *F2 just after the motions (hereinafter, referred to 
as post-motion rotation angle information where applicable) through a predetermined 



computing process based on the outputted pre-motion acceleration information and post- 
motion acceleration information; and a second computing unit for calculating the motion 
rotation angle information through a predetermined computing process based on the 
calculated pre-motion rotation angle information and post-motion rotation angle information. 
Preferably, in the first computing unit, the pre-motion rotation angle information <D1 



and the post-motion rotation angle information <D2 become <E>1 = tail 



V A bz\ 



and 



02 = tan 



-\( A byl ^ 
V A bz2 J 



, respectively, and, here, if coordinate axes of a body frame of the input 



part are denoted as X, Y, and Z, A by i and A by 2 denote the pre-motion acceleration information 
and the post-motion acceleration information for the Y axis respectively, and A bz i and A bZ 2 
denote the pre-motion acceleration information and the post-motion acceleration information 
for the Z axis respectively, and, if coordinate axes for the navigation frame are denoted as X 0 , 
Y 0 , and Zo, ¥1 and ¥2 denote the pre-motion rotation angle information and the post-motion 
rotation angle information for the Zo axis, and 01 denotes the pre-motion rotation angle 
information for a Yi axis indicating an axis after the Y 0 axis is rotated as much as ¥1, 92 
denotes the post-motion rotation angle information for Yi axis indicating an axis after the Y 0 
axis is rotated as much as *F2, Ol denotes the pre-motion rotation angle information for the 
X 2 indicating an axis after the X 0 axis is rotated as much as *F1 and 91, respectively, and 02 
denotes the pre-motion rotation angle information for the X 2 axis indicating an axis after the 
X 0 is rotated as much as ¥2 and 02, respectively. 

Preferably, in the first computing unit, the pre-motion rotation angle information 91 

f \ 

and the post-motion rotation angle information 92 become 01 = tan * 



A bxl 



K ^ A by\ 2 + A bz\ 2 J 



and 
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02 = tan 



- l \x2 



2 

bz2 J 



, respectively, and, here, if coordinate axes of a body frame of 



K 4 A by2 + A bz 

the input part are denoted as X, Y, and Z, A bx i and A bx2 denote the pre-motion acceleration 
information and the post-motion acceleration information for the X axis respectively, and 
A b yi and A by 2 denote the pre-motion acceleration information and the post-motion 
acceleration information for the Y axis respectively, and Abzi and A bz2 denote the pre-motion 
acceleration information and the post-motion acceleration information for the Z axis 
respectively, and, if coordinate axes for the navigation frame are denoted as X 0 , Y 0 , and Zo, 
*F1 and ¥2 denote the pre-motion rotation angle information and the post-motion rotation 
angle information for the Zo axis, and 01 denotes the pre-motion rotation angle information 
for a Yi axis indicating an axis after the Y 0 axis is rotated as much as *F1, and 02 denotes the 
post-motion rotation angle information for Yi axis indicating an axis after the Y 0 axis is 
rotated as much as *¥2. 

[34] Preferably, in the second computing unit, the motion rotation angle information O 

becomes <&(r ) = a * t + b , and, here, if ti denotes time just before the motions, t 2 denotes time 
just after the motions, a denotes 2 ) - <5>{t x )]/(t 2 - t x ), b denotes - a * t x + ) , and 

coordinate axes for the navigation frame are denoted as X 0 , Y 0 , and Zo, the *F denotes the 
rotation angle information for the Zo axis, the 0 denotes the rotation angle information for the 
Y! axis indicating an axis after the Y 0 axis is rotated as much as x ¥, and the O denotes the 
rotation angle information for the X 2 axis indicating an axis after the X 0 axis is rotated as 
much as and 0, respectively. 
[35] Preferably, in the second computing unit, the motion rotation angle information 0 

becomes 0{t) = c * t + d , and, here, if ti denotes time just before the motions, t 2 denotes time 
just after the motions, c denotes [&(t 2 ) - 6(t x )]/(t 2 - t x ) , d denotes - c * t, + 0{t x ) , and 
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coordinate axes are denoted as X 0 , Y 0 , and Zo, the ¥ denotes the rotation angle information 
for the Z 0 axis and the 0 denotes the rotation angle information for the Yi axis indicating an 
axis after the Y 0 axis is rotated as much as y ¥. 

[36] Further, the above object is achieved by providing the input system based on the 

three-dimensional inertial navigation system according to the present invention and having an 
input part and a host device, which detects motion position information corresponding to 
three-dimensional motions of the input part and outputs the detected motion position 
information to the host device. The input system comprises acceleration sensors for 
outputting motion acceleration information; a rotation angle information estimation- 
computing portion for estimating motion rotation angle information <D, 0, and *F based on 
acceleration information based on the gravitational acceleration separated from the outputted 
motion acceleration information; a conversion-computing unit for calculating motion position 
information based on the estimated motion rotation angle information and the outputted 
motion acceleration information; and an optimal plane-computing unit for projecting the 
motion position information onto an optimal plane. 

[37] The rotation angle information estimation-computing portion includes a separation 

unit for separating acceleration information based on the motions of the input part itself and 
acceleration information based on the gravitational acceleration from the outputted motion 
acceleration information based on a predetermined process; and a computing unit for 
calculating the motion rotation angle information through the predetermined computing 
process based on the acceleration information based on the separated gravitational 
acceleration. 

[38] Preferably, but not necessarily, the predetermined process for separating the 

acceleration information based on the gravitational acceleration from the motion acceleration 

information is to pass the motion acceleration information through a low-pass filter. 
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[39] 



Preferably, but not necessarily, in the computing unit, the motion rotation angle 



information O becomes O = tan 



-i 




, and, here, if coordinate axes of a body frame of 



the input part are denoted as X, Y, and Z, A by denotes acceleration information for the Y axis 
and A bz denotes acceleration information for the Z axis, and, if coordinate axes for a 
navigation frame are denoted as X 0 , Y 0 , and Zo, *F denotes rotation angle information for the 
Zo axis, 9 denotes rotation angle information for a Yi axis indicating an axis after the Y 0 axis 
is rotated as much as and O denotes rotation angle information for an X 2 indicating an axis 
after the X 0 is rotated as much as *F and 9, respectively. 
[40] Preferably, but not necessarily, in the computing unit, the motion rotation angle 



frame are denoted as X, Y, and Z, A bx denotes acceleration information for the X axis, A by 
denotes acceleration information for the Y axis, A bz denotes acceleration information for the 
Z axis, and coordinate axes for the navigation frame are denoted as Xo, Y 0 , and Zo, *F denotes 
rotation angle information for the Zo axis and 9 denotes rotation angle information for a Yi 
axis indicating an axis after the Y 0 is rotated as much as 
[41] Further, the above object is achieved by providing a trajectory estimation method for 

an input system based on a three-dimensional inertial navigation system and having an input 
part and a host device according to the present invention, which detects motion position 
information corresponding to three-dimensional motions of the input part and outputs the 
detected motion position information to the host device. The method comprises steps of (a) 
outputting motion acceleration information, pre-motion acceleration information, and post- 
motion acceleration information just after the motions; (b) estimating motion rotation angle 
information G>, 9, and ¥ through a predetermined computing process based on the outputted 



information 9 becomes 6 = tail 



-i 




and, here, if coordinate axes for the body 



pre-motion acceleration information and post-motion acceleration information; (c) calculating 
the motion position information based on the estimated motion rotation angle information and 
the outputted motion acceleration information; and (d) projecting the motion position 
information onto an optimal plane. 

[42] The step (b) includes steps of (bl) calculating pre-motion rotation angle information 

Ol, 91, and and post-motion rotation angle information 02, 02, and *F2 through a 
predetermined computing process based on the outputted pre-motion acceleration information 
and post-motion acceleration information; and (b2) calculating the motion rotation angle 
information through a predetermined computing process based on the calculated pre-motion 
rotation angle information and post-motion rotation angle information. 

[43] Preferably, but not necessarily, in the step(bl), the pre-motion rotation angle 

information Ol and the post-motion rotation angle information 02 become 



<Di = tan" 



and 02 = tan 



f A 

^by2 



respectively, and, here, if coordinate axes of a 



V ) 

body frame of the input part are denoted as X, Y, and Z, A by i and A by 2 denote the pre-motion 
acceleration information and the post-motion acceleration information for the Y axis 
respectively, and A bz i and A bZ 2 denote the pre-motion acceleration information and the post- 
motion acceleration information for the Z axis respectively, and, if coordinate axes for the 
navigation frame are denoted as X 0 , Y 0 , and Zo, *F1 and ¥2 denote the pre-motion rotation 
angle information and the post-motion rotation angle information for the Z 0 axis, and 91 
denotes the pre-motion rotation angle information for a Yi axis indicating an axis after the Y 0 
axis is rotated as much as ¥1, 92 denotes the post-motion rotation angle information for Yi 
axis indicating an axis after the Y 0 axis is rotated as much as ¥2, Ol denotes the pre-motion 
rotation angle information for the X 2 indicating an axis after the X 0 axis is rotated as much as 
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¥1 and 81, respectively, and 02 denotes the pre-motion rotation angle information for the X 2 
axis indicating an axis after the X 0 is rotated as much as *F2 and 02, respectively. 

Preferably, but not necessarily, in the step(bl), the pre-motion rotation angle 
information 91 and the post-motion rotation angle information 02 become 



« = tan" 1 



f \ ( n 



\x\ 



and 02 = tan 



/— f— — toil , 2 



\x2 



2 

*b Z 2 j 



, respectively, and, here, if 



coordinate axes of a body frame of the input part are denoted as X, Y, and Z, A bx i and A bX 2 
denote the pre-motion acceleration information and the post-motion acceleration information 
for the X axis respectively, and A by i and A by 2 denote the pre-motion acceleration information 
and the post-motion acceleration information for the Y axis respectively, and A bz i and A bz2 
denote the pre-motion acceleration information and the post-motion acceleration information 
for the Z axis respectively, and, if coordinate axes for the navigation frame are denoted as X 0 , 
Y 0 , and Zo, *F1 and *F2 denote the pre-motion rotation angle information and the post-motion 
rotation angle information for the Zo axis, and 01 denotes the pre-motion rotation angle 
information for a Yi axis indicating an axis after the Yo axis is rotated as much as ¥1, and 02 
denotes the post-motion rotation angle information for Yi axis indicating an axis after the Y 0 
axis is rotated as much as *F2. 

Preferably, but not necessarily, in the step(b2), the motion rotation angle information 
O becomes 0(r ) = a*t + b and, here, if ti denotes time just before the motions, t 2 denotes 
time just after the motions, a denotes [<I>(f 2 ) - )]/(r 2 - t x ), b denotes - a * t , + ®{t x ), and 

coordinate axes are denoted as X 0 , Y 0 , and Zo, the ¥ denotes the rotation angle information 
for the Zo axis, the 0 denotes the rotation angle information for the Yi axis indicating an axis 
after the Y 0 axis is rotated as much as V F, and the O denotes the rotation angle information for 
the X 2 axis indicating an axis after the X 0 axis is rotated as much as *F and 0, respectively. 
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[46] Preferably, but not necessarily, in the step(b2), the motion rotation angle information 

0 becomes 0{t ) = c*t + d and, here, if ti denotes time just before the motions, t2 denotes time 
just after the motions, c denotes [0(t 2 ) - 0(t x )]/(t 2 - t l ) , d denotes —c*t t + 0{t x ) , and 

coordinate axes are denoted as X 0 , Y 0 , and Zo, the *F denotes the rotation angle information 
for the Zo axis and the 0 denotes the rotation angle information for the Yi axis indicating an 
axis after the Yo axis is rotated as much as V P. 
[47] Further, the above object is achieved by providing a trajectory estimation method for 

an input system based on a three-dimensional inertial navigation system and having an input 
part and a host device, which detects motion position information corresponding to motions 
of the input part based on three-dimensional motions of the input part and outputs the 
detected motion position information to the host device. The method comprises steps of: (a) 
outputting motion acceleration information; (b) estimating motion rotation angle information 
O, 0, and *F based on acceleration information based on the gravitational acceleration 
separated from the outputted motion acceleration information; (c) calculating motion position 
information based on the estimated motion rotation angle information and the outputted 
motion acceleration information; and (d) projecting the motion position information onto an 
optimal plane. 

[48] The step (b) includes steps of (bl) separating acceleration information based on the 

motions of the input part itself and acceleration information based on the gravitational 
acceleration from the outputted motion acceleration information based on a predetermined 
process; and (b2) calculating the motion rotation angle information through the 
predetermined computing process based on the acceleration information based on the 
separated gravitational acceleration. 
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[49] 



Preferably, but not necessarily, the predetermined process in the step (bl) is to pass 



the motion acceleration information through a low-pass filter. 



[50] 



Preferably, but not necessarily, the motion rotation angle information O in the step 



(b2) becomes <t> = tail 



-i 




and, here, if coordinate axes of a body frame of the input 



part are denoted as X, Y, and Z, A by denotes acceleration information for the Y axis and A bz 
denotes acceleration information for the Z axis, and, if coordinate axes for a navigation frame 
are denoted as Xo, Yo, and Z 0 , denotes rotation angle information for the Zo axis, 0 denotes 
rotation angle information for a Yi axis indicating an axis after the Y 0 axis is rotated as much 
as and O denotes rotation angle information for an X 2 indicating an axis after the X 0 is 
rotated as much as *F and 0, respectively. 
[51] Preferably, but not necessarily, the motion rotation angle information 0 in the step 



denoted as X, Y, and Z, A bx denotes acceleration information for the X axis, A by denotes 
acceleration information for the Y axis, A bz denotes acceleration information for the Z axis, 
and coordinate axes for the navigation frame are denoted as Xo, Y 0 , and Zo, ¥ denotes 
rotation angle information for the Zo axis and 0 denotes rotation angle information for a Yi 
axis indicating an axis after the Y 0 is as much as x ¥. 
[52] Accordingly, the present invention can recover the three-dimensional motions by 

using only the three acceleration sensors without the angular velocity sensors in the three- 
dimensional inertial navigation system. 



(b2) becomes 0 = \fxx\ 



and, here, if coordinate axes for the body frame are 




BRIEF DESCRIPTION OF THE DRAWINGS 



[53] The features and advantages of the present invention will become more apparent by 

describing in detail exemplary embodiments thereof with reference to the following 
drawings in which like reference numerals refer to like elements, and wherein: 

[54] Fig. 1 is a block diagram for showing an input system using a conventional inertial 

navigation system; 

[55] Fig. 2 is a block diagram for showing an input system including an input part 

equipped with a rotation angle information estimation-computing portion according to an . 

embodiment of the present invention; 
[56] Fig. 3 is a block diagram for showing an input system including an input part 

equipped with a rotation angle information estimation-computing portion according to 

another embodiment of the present invention; 
[57] Fig. 4A is a graph for showing acceleration information on motions outputted from 

acceleration sensors; 

[58] Fig. 4B is a graph for showing acceleration information in case that only the 

acceleration information based on the gravitational acceleration is separated from the 
acceleration information of Fig. 4A when a low-pass filter is provided in a separation unit; 

[59] Fig. 5 is a flow chart for showing a trajectory estimation method for an input system 

of Fig. 2 according to an embodiment of the present invention; 

[60] Fig. 6 is a flow chart for showing a trajectory estimation method for an input system 

of Fig. 3 according to another embodiment of the present invention; 

[61] Fig. 7 A is a view for showing results of certain numbers written by an input system to 

which the present invention is applied; 

[62] Fig. 7B is a view for showing results of other certain numbers written by an input 

system to which the present invention is applied; and 
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Fig. 7C is a view for showing results of certain letters written by an input system to 
which the present invention is applied. 

DETAILED DESCRIPTION OF ILLUSTRATIVE, NON-LIMITING EMBODIMENTS OF 

THE INVENTION 

Hereinafter, an input system according to an exemplary embodiment of the present 
invention is described in detail with reference to the accompanying drawings. 

Fig. 2 is a block diagram for showing an input system including an input part 100 
equipped with a rotation angle information estimation -computing portion 120 according to an 
embodiment of the present invention. 

The input system according to the present invention has the input part 100 and a host 
device 300. 

The input part 100 includes acceleration sensors 110, the rotation angle information 
estimation-computing portion 120, a conversion-computing unit 130, an optimal plane- 
computing unit 135 and a transmission unit 140. 

The acceleration sensors 1 10 are provided for three X, Y, and Z axes of a body frame 
of the input part 100 according to three x, y, and z axes to be based on motions of the input 
part 100. The acceleration sensors 1 10 detect and output acceleration information on motions 
(hereinafter, referred to as motion acceleration information where applicable), acceleration 
information just before the motions (hereinafter, referred to as pre-motion acceleration 
information where applicable), and acceleration information just after the motions 
(hereinafter, referred to as post-motion acceleration information where applicable), 
respectively. It is preferable for the three acceleration sensors to be provided in the input part 
100, but it is possible to provide two acceleration sensors that can detect acceleration 
information on two axes. 
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[69] Definitions for the acceleration information on motions, acceleration information just 

before the motions, and acceleration information just after the motions, are as follows. 

[70] In order to do trajectory estimations for the input part 100 according to an 

embodiment of Fig. 2, an assumption is necessary that the input part 100 has to remain still 
just before and after intended motions to be taken by the input part 100. Accordingly, the 
acceleration sensors 110 according to the present invention can detect acceleration 
information just before and after intended motions to be taken by the input part 100. The pre- 
motion acceleration information refers to acceleration information just before intended 
motions. Further, the post-motion acceleration information refers to acceleration information 
just after intended motions. Motion acceleration information refers to acceleration 
information based on intended motions to be taken by a user. 

[71] The rotation angle information estimation-computing portion 120 computes and 

estimates rotation angle information based on pre-motion acceleration information and post- 
motion acceleration information outputted from the acceleration sensors 110. 

[72] The rotation angle information estimation-computing portion 120 according to an 

embodiment of the present invention includes a first computing unit 121 and a second 
computing unit 122. 

[73] The first computing-unit 121 inputs pre-motion acceleration information and post- 

motion acceleration information on the input part 100 from the acceleration sensors 110. 

[74] The first computing unit 121 computes O and 0 of rotation angle information just 

before motions through a predetermined computing process based on pre-motion acceleration 
information. Here, the rotation angle information just before motions is pre-motion rotation 
angle information corresponding to the pre-motion acceleration information. 

[75] The first computing unit 121 computes <I> and 0 of rotation angle information just 

after motions through a predetermined computing process based on post-motion acceleration 
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information. In here, the rotation angle information just after motions is post-motion rotation 
angle information corresponding to the post-motion acceleration information. 

In case that coordinate axes for the body frame are denoted as X, Y, and Z, 
acceleration information for the X axis of the body frame is denoted as A bx , acceleration 

information for the Y axis of the body frame as A by , acceleration information for the Z axis 

of the body frame as A bz , rotation angle information for a Z 0 axis as \|/, and rotation angle 

information for a Y x axis, indicating an axis after a Y 0 axis is rotated as much as *F, as 0, 

Equation 4 expresses rotation angle information O for a X 2 axis indicating an axis after the 

X 0 axis is rotated as much as *F and 0, respectively. 

[Equation 4] 



In case that coordinate axes for the body frame are denoted as X, Y, and Z, 
acceleration information for the X axis of the body frame is denoted as A bx , acceleration 

information for the Y axis of the body frame as A by , acceleration information for the Z axis 

of the body frame as A bz , and rotation angle information for a Z 0 axis as V P, rotation angle 

information 0 for Y l axis indicating an axis after a Y 0 axis is rotated as much as *F can be 

expressed in Equation 5 as follows: 

[Equation 5] 



Equations 4 and 5 are formulae by which the <J> and 0 of rotation angle information 
can be respectively calculated from acceleration information at the state of no motion. 
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[79] The second computing unit 122 inputs the O and 9 of pre-motion rotation angle 

information computed from the first computing unit 12 1. 
[80] The second computing unit 122 inputs the O and 9 of post-motion rotation angle 

information computed from the first computing unit 121. 
[81] The second computing unit 122 computes the motion rotation angle information O 

through a predetermined computing process based on the O of the inputted pre-motion 

rotation angle information and the Oof the inputted post-motion rotation angle information. 
[82] The second computing unit 122 computes the motion rotation angle information 9 

through a predetermined computing process based on the 9 of the inputted pre-motion 

rotation angle information and the 9 of the inputted post-motion rotation angle information. 
[83] If the time just before a motion is denoted as t x , the time just after a motion as t 2 , 

[0>(r 2 ) - <b(t x )]/(t 2 - t x ) as a, and - a * t l + Ofo ) as b, the O(r) of the motion rotation angle 

information can be expressed in Equation 6 as follows: 

[Equation 6] 

0(r) = a*f + b 

[84] Further, if the time just before a motion is denoted as t x , the time just after a motion 

as t 2 , [6(t 2 ) - 0{t x )]/(t 2 -t x ) as c, and - c * t x + &(t x ) as d, the 0(t) of the motion rotation 
angle information can be expressed in Equation 7 as follows: 

[Equation 7] 

6{t) = c*t + d 

[85] The conversion-computing unit 130 inputs motion acceleration information from the 

acceleration sensors 1 10, and inputs motion rotation angle information estimated from the 
rotation angle information estimating-computing portion 120. Based on the inputted motion 
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acceleration information and motion rotation angle information, motion velocity information 

V n and motion position information P n for the navigation frame are calculated as in the 

conversion-computing unit of Fig. 1. 
[86] The optimal plane-computing unit 135 projects onto a virtual two-dimensional 

optimal plane the motion position information outputted from the conversion-computing unit 

130 to obtain coordinate values. 
[87] The transmission unit 140 outputs to the host device 300 the projected position 

information outputted from the optimal plane-computing unit 135. The transmission unit 140 

can be constructed with a wired communication module, but, preferably, with a wireless 

communication module such as Bluetooth®. 
[88] The host device 300 displays position information corresponding to motions based on 

the motions of the input part 100. The host device 300 is preferably a personal portable 

device. 

[89] Fig. 3 is a block diagram for showing an input system including the input part 200 

equipped with the rotation angle information estimation-computing unit 120 according to 
another embodiment of the present invention. 

[90] As described in Fig. 2, the input system according to the present invention includes 

the host device 300 and the input part 200. 

[91] The input part 200 has acceleration sensors 210, a rotation angle information 

estimation-computing portion 220, a conversion-computing unit 230, an optimal plane- 
computing unit 235, and a transmission unit 240. 

[92] The description on the conversion-computing unit 230 and the transmission unit 240 

is the same as in Fig. 2. 

[93] Accordingly, hereinafter, a description is made only of the acceleration sensors 210 

and the rotation angle information estimation-computing portion 220. 
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[94] The acceleration sensors 210 according to the present embodiment detect and output 

motion acceleration information. Differently from the embodiment of Fig. 2, there is no need 
to detect acceleration information just before and after motions. 

[95] The rotation angle information estimation-computing portion 220 according to the 

present embodiment includes a separation unit 221 and a computing unit 222. 

[96] The separation unit 221 inputs motion acceleration information outputted. Thus, the 

separation unit 221 separates acceleration information based on the motion of the input part 
200 itself and acceleration information based on the gravitational acceleration from the 
inputted motion acceleration information through a predetermined process. 

[97] For the predetermined process, it is preferable to provide a low-pass filter (LPF) to the 

separation unit 221. 

[98] In general, the acceleration information based on the gravitational acceleration exists 

in a lower frequency bandwidth than the acceleration information based on motions 
themselves. Accordingly, in case that a low-pass filter is provided to the separation unit 221, 
the acceleration information based on the gravitational acceleration is filtered by the 
separation unit 221. 

[99] The computing unit 222 is inputted with acceleration information based on the 

gravitational acceleration. 
[100] The computing unit 222 calculates motion rotation angle information through 

Equations 4 and 5 based on acceleration information based on the gravitational acceleration. 
[101] In general, a stationary object has no motions and is affected only by the gravitational 

force, so that acceleration information based on the gravitational acceleration of the motion 

acceleration information corresponds to a stationary state. 

[102] Further, as described above, Equations 4 and 5 can be applied only in a stationary 

state. Therefore, Equations 4 and 5 can be applied to the acceleration information based on 
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the gravitational acceleration. The computing unit 222 inputs the acceleration information 
based on the gravitational acceleration, and calculates motion rotation angle information 
through applications of Equations 4 and 5. 

[103] The calculated motion rotation angle information and the motion acceleration 

information outputted from the acceleration sensors 210 are transmitted to the conversion- 
computing unit 230, and the motion position information for the navigation frame calculated 
in the conversion-computing unit 230 is transmitted to the host device 300 through the 
optimal plane-computing unit 135. 

[104] However, the present disclosure does not perform any estimation for \\f. This is 

because the \|/ itself does not affect the trajectory recovery for the input part. Therefore, the 
values for the \\f in the two embodiments are set to zero. 

[105] Fig. 4 A is a graph for showing motion acceleration information outputted from the 

acceleration sensors 1 10. 

[106] In general, motion acceleration information outputted from the acceleration sensors is 

composite information of acceleration information based on the gravitational acceleration and 
acceleration information based on motions themselves. As shown in Fig. 4A, acceleration 
information based on motions themselves generally exists in a frequency band ranging from 
1Hz to 20Hz, and acceleration information based on the gravitational acceleration exists in a 
frequency band near 0Hz. 

[107] Fig. 4B is a graph for showing acceleration information in case that only the 

acceleration information based on the gravitation acceleration is separated from the 
acceleration information of Fig. 4A with the separation unit 221 equipped with a low-pass 
filter. 

[108] In case that the separation unit 221 has a low-pass filter, acceleration information in a 

high-frequency band is not filtered for an output, but only the acceleration information in a 
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low-frequency band is outputted. Therefore, as shown in Fig. 4B, only the acceleration 

information based on the gravitational acceleration is outputted. 
[109] Fig. 5 is a flow chart for showing a trajectory estimation method based on the input 

system according to the embodiment of Fig. 2. 
[110] Descriptions are made of the method with reference to Fig. 2 and Fig. 5. First, a user 

takes certain motions with the input part 100. The acceleration sensors 1 10 in the input part 

100 detect and output motion acceleration information based on the motions of the input part 

100 (S400). 

[Ill] The outputted acceleration information consists of pre-motion acceleration 

information, motion acceleration information, and post-motion acceleration information. The 
outputted pre-motion and post-motion acceleration information is inputted to the rotation 
angle information estimation-computing portion 120. The first computing unit 121 in the 
rotation angle information estimation-computing portion 120 applies Equations 4 and 5 to the 
pre-motion acceleration information and the post-motion acceleration information. Thus, 
pre-motion rotation angle information corresponding to the pre-motion acceleration 
information is calculated, and post-motion rotation angle information corresponding to the 
post-motion acceleration information is calculated (S410). 

[112] The calculated pre-motion rotation angle information and the post-motion rotation 

angle information are outputted to the second computing unit 122. The second computing 
unit 122 applies Equations 6 and 7 to the pre-motion rotation angle information and the post- 
motion rotation angle information. Accordingly, the motion rotation angle information is 
calculated (S420). 

[113] The calculated motion rotation angle information is outputted to the conversion- 

computing unit 130. Further, the acceleration information outputted from the acceleration 

sensors 110 is outputted to the conversion-computing unit 130. The conversion-computing 
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unit 130 calculates motion position information for the navigation frame based on the 
inputted acceleration information and angular velocity information (S430). 
[1 14] The motion position information for the calculated navigation frame is projected onto 

an optimal plane by the optimal plane-computing unit 135 (S440). The position information 
projected on the optimal plane is outputted to the host device 300 through the transmission 
unit 140 (S450). 

[115] Fig. 6 is a flow chart for showing a trajectory estimation method for the input system 

according to the embodiment of Fig. 3. 

[116] Descriptions are made with reference to Fig. 3 and Fig. 6. First, a user takes certain 

motions with the input part 200. The acceleration sensors 210 in the input part 200 outputs 
motion acceleration information based on the motions of the input part 200 (S500). 

[1 17] The outputted motion acceleration information consists of acceleration information 

based on the motions themselves and acceleration information based on the gravitational 
acceleration. The outputted motion acceleration information is inputted to the rotation angle 
information estimation-computing portion 220. The separation unit 221 in the rotation angle 
information estimation-computing portion 220 separates the inputted motion acceleration 
information into the acceleration information based on the motions themselves and the 
acceleration information based on the gravitational acceleration (S510). 

[118] The acceleration information based on the separated gravitational acceleration is 

inputted to the rotation angle information estimation-computing portion 220. The computing 
unit 222 in the rotation angle information estimation-computing portion 220 applies 
Equations 4 and 5 to the acceleration information based on the gravitational acceleration of 
the separated acceleration information. Accordingly, motion rotation angle information is 
calculated (S520). 
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[119] The calculated motion rotation angle information is outputted to the conversion- 

computing unit 230. Further, the acceleration information outputted from the acceleration 
sensors 210 is outputted to the conversion-computing unit 230. The conversion-computing 
unit 230 calculates motion position information for the navigation frame based on the 
inputted acceleration information and angular velocity information (S530). 

[120] The calculated motion position information for the gravitational frame is projected 

onto an optimal plane by the optimal plane-computing unit 235 (S540). The position 
information projected on the optimal plane is outputted to the host device 300 through the 
transmission unit 240 (S550). 

[121] Fig. 7 A illustrates results of certain numbers written by an input apparatus to which 

the present invention is applied. 

[122] Fig. 7 A shows the experimental results for numbers 0-3. 

[123] The one-dotted lines in Fig. 7 A indicate numbers written by an input part without 

angular velocity sensors according to the embodiment of Fig. 2. The dotted lines in Fig. 7A 
indicate numbers written by the input part without angular velocity sensors according to the 
embodiment of Fig. 3. The solid lines indicate numbers directly written on a tablet by a 
stylus pen. The two-dotted lines in Fig. 7A indicate numbers written by an input part with 
angular velocity sensors. 

[124] Fig. 7B illustrates results of other certain numbers written by the input part to which 

the present invention is applied, and Fig. 7C illustrates results of certain letters written by the 

input part to which the present invention is applied. 
[125] Fig. 7B shows an experimental result for numbers 6-9, and Fig. 7C shows an 

experimental result for English letters G, H, I, and J. 
[126] The descriptions on the one-dotted lines, dotted lines, solid lines, and two-dotted lines 

are the same as in Fig. 7A. 
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[127] As shown in Fig. 7A, Fig. 7B, and Fig. 7C, the numbers and letters that are written by 

the input part with angular velocity sensors and directly written on the tablet are nearly the 
same as those that are written by the input parts according to the embodiments of the present 
invention, and, further, it can be seen that the degrees to which the writings can be recognized 
are nearly the same 

[128] However, the embodiment presented in Fig. 2 and the embodiment present in Fig. 3 

are each structured to perform a predetermined computing process in the input part so that 
only the position information as final processing result values is transmitted to the host device. 

[129] However, the above structure is exemplified for the above embodiments, and, for a 

different embodiment, the structure may be constructed in a different way in order that the 
input part detects only acceleration information outputted from the acceleration sensors, the 
detected acceleration information is transmitted to the host device through the transmission 
unit, and the host device performs the computing process by the rotation angle information 
estimation-computing portion, the computing process by the conversion-computing unit, and 
the optimal plane-computing process. 

[130] Further, for another embodiment, the structure may be constructed so that the input 

part detects acceleration information outputted from the acceleration sensors, the rotation 
angle information estimation-computing portion calculates rotation angle information based 
on the detected acceleration information, the calculated rotation angle information is 
transmitted to the host device through the transmission unit, and the host device performs the 
computing process by the conversion-computing unit and the optimal plane-computing 
process. 

[131] The input system based on the three-dimensional inertial navigation system and a 

trajectory estimation method thereof according to the present invention can recover three- 
dimensional motions in use of only three acceleration sensors without angular velocity 
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sensors, rather than using three acceleration sensors and three angular velocity sensors as in 
the three-dimensional inertial navigation system. Further, the present invention can input any 
writings a user desires without limit to a time and place. Accordingly, the input system can 
be manufactured at a low cost since only three acceleration sensors are used, manufactured 
light-weight to be convenient in carrying, reduce power consumption, and eliminate the 
troublesomeness of initial correction work. 
[132] While the invention has been shown and described with reference to exemplary 

embodiments thereof for illustrative purposes, it will be understood by those skilled in the art 
that various changes in form and details may be made therein without departing from the 
spirit and scope of the invention as defined by the appended claims. 
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